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ABSTRACT 

We discover a pronounced dependence of the strength of the soft X-ray orbital modulation 
and the spectral hardness in Cyg X-l in the hard state on its superorbital phase. We find our 
results can be well modelled as a combination of two effects: the precession of the accretion 
disc (which appears to cause the superorbital flux modulation) and the orbital-phase depen- 
dent X-ray absorption in an accretion bulge, located at the accretion disc edge close to the 
supergiant companion but displaced from the line connecting the stars by about 25°. Our 
findings are supported by the distribution of the X-ray dips showing concentration towards 
zero superorbital phase, which corresponds to the bulge passing through the line of sight. We 
Fourier analyse our model, and find it explains the previous finding of asymmetric beat (be- 
tween the orbital and superorbital modulations) frequencies in the observed power spectrum, 
provided the disc precession is prograde. On the other hand, we find no statistically significant 
changes of the orbital modulation with the superorbital phase in the 15-GHz radio data. This 
absence is consistent with the radio being emitted by a jet in the system, in which case the 
orbital modulation is caused by wind absorption far away from the disc. We also find that both 
the X-ray and radio fluxes of Cyg X-l in the hard state on time scales >10 4 -s have lognormal 
distributions, which complements a previous finding of a lognormal flux distribution in the 
hard state on ^1-s time scales. We point out that the lognormal character of the flux distribu- 
tion requires that flux logarithms rather than fluxes themselves should be used for averaging 
and error analysis. We also provide a correct formula for the uncertainty of rms of a light 
curve for the case when the uncertainty is higher than the measurement. 

Key words: accretion, accretion discs - radio continuum: stars - stars: individual: Cyg X-l 
- stars: individual: HDE 226868 - X-rays: binaries - X-rays: stars. 



1 INTRODUCTION 

A number of X-ray binaries show flux periodicities at their respec- 
tive orbital period, which may be caused by a number of effects. 
First, the source associated with the compact object in a binary may 
be eclipsed by th e companion (usually of high mass) (see e.g. a list 
I Wen et al. 2006). Second, a flux modulation may be caused by an 
optically-thick disc rim (which is highest at the point of impact of 
the gas stream from the inner Langrangial point in case of a donor 
filling its Roche lobe) , obscuring the disc and/ or its corona (e.g., 
IWhite & Swanklll999l : iHellier & Masor]|l989l) . This obscuration 
may lead to strong partial eclipses in so-called X-ray dippers. More 
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generally, the disc and any associated structures may depart from 
its axial symmetry due to the influence of the companion, which 
may cause an orbital modulation. Third, wind from a high-mass 
companion may absorb/scatter the emission from the vicinity of 
the compact object, and the degree of absorption will depend on the 
orbital phase. In the case of Cyg X-l, both X-ray and radio emis- 
sion are mo dulated by this ef fect, which modulations were m od- 
elled by, e.g. JWen et aj]jl999h and lSzostek & Zdziarskild2007l) . re- 
spectively. Fourth, phase-dependent absorption (via photon-photon 
pair production) of high-energy 7-rays may occur in a photon field 
axially asymmetric with respect to the compact object, especially 
that of the stellar photons (e.g., lBedna rek 2006). A fifth effect of 
the companion is reflection or reprocessing of the emission from 
around the compact object on the surface of the companion fac- 
ing the compact object. This effects appears to be responsible for, 
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e.g., the UV flux modulation from the X-ray binary 4U 1820- 
303 dArons & Kin3ll993l : lAnderson etal] [l997). Finally, the op- 
tical/UV emission of the companion will be modulated if its shape 
departs from the spherical symmetry by partially or fully filling its 
Roche lobe, which effect is seen in Cyg X-l, e.g.. lBrocksopp et at] 
dl999bl) . 

Then, there will be an intrinsic dependence of the emitted flux 
on the orbital phase if the orbit is elliptical. This leads, e.g., to pe- 
riodic outbursts at the periastro n of Cir X-l jParkinson et alj2003h 
and Be/X-ray binaries (see, e.g., Coe 2000; Negueruela 2004 for re- 
views) in X-rays, and sometimes, at other wavelengths. Also, some 
orbital flux modulation may be due to the Dopple r effect, which is 
in principle observable jPostnov & Sha kura 1987), but has not yet 
been detected in a binary. (Obviously, the Doppler effect leads to 
widely observed shifts of spectral lines from binaries.) 

In addition, a number of X-ray binaries show modulation 
at periods much longer than their orbital peri ods, so-ca l led su - 
perorbital periodicity, see, e.g., a partial list in IWen et al. I d2006h . 
In particular, Cy g X-l shows such periodicity with the period 
of ~150 d (e.g., lBrocksoppetal.lll999j:lKaritskavaetal.ll200ll ; 
Ozdemir & Demircan 2001; Lachowicz et al. 2006, hereafter L06; 
Ibragimov. Zdziarski & Poutanenll 20071 . hereafter Paper I). The ob- 
served superorbital variability appears in most cases compatible 
with being caused by accretion disc and/or jet precession, which 
either results in variable obscuration of emitted X-rays as in Her 
X-l dKatzll 19731) . or changes th e viewing angle of the presumed 
anisotropic emitter, as in SS 433 dKatzll 980) or Cyg X-l (e.g., L06, 
Paper I), or both. The only known exception, in which the super- 
orbital periodicity is clearly caused by modulation of the accretion 
rate (and t hus not by a changing viewing angle of the source), is 4U 
1820-303 dZdziarski. Wen & Gierlinskil2007ah . 

A number of binaries show both orbital and superorbital mod- 
ulations. Those currently known are LMC X-4, 2S 01 14+650, SMC 
X-l, Her X-l, SS 433, 4U 1820-303 and Cyg X-l. An interesting 
issue then is whether there is any dependence of the parameters 
of the orbital modulation on the superorbital phase (or, similarly, 
on an average of the flux level). The shape of the profile of the 
orbital mod ulation in Her X-l wa s found to depend on its superor- 
bital phase dScott & Leahvi ri999). which appears to be due to the 
shadowing effect of the precessing accretion disc and scattering in 
its wind in that system. Recently, analogous dependencies of the 
shape of the orbital modulation on the average flux level have been 
found in LMC X-4, S MC X- l, Her X-l, as well as in Cen X-3 
( Raic hur & Pauill2008h . Then, IZdziarski et all [2007b) found such 
a dependence in 4U 1820-303 (of both the amplitude and the phase 
of the minimum flux) and interpreted it in terms of the size of the 
disc rim (partially obscuring the central source) changing with the 
variable accretion rate. 

In addition, there is the case of the peculiar Be/X-ray bi- 
nary LS I +61° 303, which shows orbital variability in the ra- 
dio, X-ray and TeV emission, and a superorbital variability of the 
peak radio flux during an orbit ( Grego ry? Peracaula & Taylor! 19991 ; 
lGregorvll2002l) . lGregorvl d2002h found a marked dependence of 
the phase of the peak of the orbital radio modulation on the su- 
perorbital phase in LS I +61°303. The presence of such a de- 
pendence may be due to interaction of the pulsa r in that system 
with a variable circumstellar Be decretion disc dGregorvl 12002; 
IZdziarski. Neronov & Chern vakova 20081). 

It is of considerable interest to find out whether orbital modu- 
lation depends on the superorbital phase in Cyg X-l, the archetypi- 
cal and very well studied black-hole syste m with a high-m ass com- 
panion, the OB supergiant HDE 226868 dWalbornlll973l) . In this 
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Figure 1. (a) The histogram of the count rates (dwell-by-dwell data) ob- 
served from Cyg X-l for the ASM A, B, and C detector channels, (b) Dis- 
tribution of the hardness ratios B/A, C/B, and C/A. (c) Distribution of the 
flux in the Ryle data. The flux units for the ASM and Ryle data are count 
s _1 and Jy, respectively. The solid curves give the best- fitting lognormal 
distributions. 

work, we study this issue and find that such dependence exists and 
is very strong in soft X-rays. We then explain it theoretically in 
terms of orbital-phase dependent absorption in the stellar wind in- 
teracting with the outer accretion disc. 



2 THE LIGHT CURVES AND THEIR ANALYSIS 
2.1 Data 

We use the X-ray dwell data (MJD 50087-53789, i.e., 1996 Jan- 
uary 5-2006 February 23; note a misprint in the start date in Pa- 
per I) obtained with three Scanning Shadow Cameras of the All- 
Sky Monitor (ASM) aboard Rossi X-ray Timing Explo rer (RXTE; 
iBradt. Rothschild & Swaiikl 1 19931 ; iLevine et ali Il996h . with the 
channels A, B, and C corresponding to the photon energy inter- 
vals of 1.5-3 keV, 3-5 keV, and 5-12 keV, respectively. We also 
use the corresponding 15-GHz radio data from the Ryle Tele- 
scope of the Mullard Radio As tronomy Observatory (see, e.g., 
Poole v. Fender & Brocksopd[l999l ; L06 for earlier analyses of the 
observations of Cyg X-l). 

Because Cyg X-l is a highly variable source and the effects 
we search for are rather weak, we need to select accurately a ho- 
mogeneous set of data. For most of the analysis in the paper, we 
use the data corresponding to the hard spectral state following the 
criteria defined in section 2 of Paper I. We require the average pho- 
ton spectral index derived from the RXTE/ 'ASM fluxes to be < 2.1 
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( Zdzi arski et al.| [2002). and additionally we exclude hard-state in- 
tervals with high X-ray variability, namely we include only those 
30-d intervals of the ASM data where < 40 per cent of points ex- 
ceed by 4<t the average flux in the reference interval MJD 50660- 
50990. This has resulted in considering the following time in- 
tervals: MJD 50350-50590, 50660-50995, 51025-51400, 51640- 
51840, 51960-52100, 52565-52770, 52880-52975, 53115-53174, 
53554-53690 (see fig. 1 in Paper I). 



2.2 Mean fluxes and variance 

We generality follow the method of analyzing light curves de- 
scribed in Paper I, but with some modifications necessitated by the 
sc ientific goal of the prese nt work. We use the orbital ephemeris 
of lBrocksopp et alj {l999b) and the superorbital ephemeris of L06, 
see equations (1) and (4), respectively, in Paper I. We use the val- 
ues of the orbital and superorbital periods of P = 5.599829 d 
and P sup = 151.43 d, respectively. We first divide an analyzed 
light curve into bins with the length of P/20. Then we average all 
points falling into a given bin weighted by the inverse squares of 
their measurement errors, obtaining the binned light curve, Fi. In 
this way, we avoid any contribution to our folded/averaged light 
curves from the source variability on time scales shorter than that 
corresponding to the length of our chosen phase bin (see Paper I). 
Note that unlike the method in Paper I, we do not prewhiten the 
light curves, i.e., do not subtract variability at one period in order 
to detect more clearly variability at another period. 

We have then looked into statistical properties of our distribu- 
tions. We plot histograms of the fluxes for the ASM and Ryle data 
in Fig. [TJ We see that each of the histograms follows a lognormal 
distribution and it is completely inconsistent with a normal one. 
Our finding of the lognormal form of the variability of Cyg X-l in 
the hard state on long time scales ( ^1/10-d to yr) in both X-rays 
and ra dio is supplemental to that of lUttlev. McHardv & Vaughanl 
(2005), who found the same type of distribution in X-rays on short 
time scales, ~0.1-10 s, also in the hard state. This form of the flux 
distribution has important implications for calculating flux aver- 
ages and the intrinsic dispersion, i.e., the standard deviation in the 
data. Namely, the standard-deviation error estimate based on the 
rms, namely 
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2 E;=i 0e» . 
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provides an unbiased estimate of the true standard-deviation er- 
ror of the average of xt only if the distribution of Xi is normal 
( iBevington & Robinson! fl992i) . Therefore, for the purpose of cal- 
culating the averages and the rms standard deviations for our light 
curves, we have converted the count rates or fluxes in our binned 
light curves, Fi, into its logarithm, Gi = In Fi, with Gi having 
now the distributions close to normal. 

We then separate the light curves binned based on the orbital 
phase into superobital phase bins of the length of P sup /8, with the 
mid-point of the first and the fifth bin at $ = and 0.5, respec- 
tively. Here, either the orbital phase, <f>, or superorbital phase, is 
defined in the 0-1 interval, and corresponds to the flux minimum 
as defined by the respective ephemeris. Then, we calculate folded 
and averaged profiles (of Gi = \nFi) of the orbital modulation 
within each superorbital phase bin, i.e., 
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where i £ (j, k) counts over all points, i, falling into a given su- 
perorbital bin, j, and the orbital bin, k, and Ijk is the number of 
such points. We estimate the error of this average using equation 
{Q, i.e., 
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Note that this error estimate accounts for both the aperiodic vari- 
ability of the source, i.e., intrinsic dispersion of individual fluxes 
contributing to a given orbital/superorbital bin (usually dominat- 
ing), and the dispersion due to measurement errors. Also, since we 
use logarithms, ajk represents a fractional error (and should not be 
divided by Gjk). 

The average and the average square error in a given superor- 
bital bin are, 
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respectively, where K = 20 is the number of orbital bins. 

We need to characterize the strength of a given modulation. 
One way of doing it without making any assumptions about its 
shape is to measure the fractional rms of a given orbital modula- 
tion profile. To do it, we calculate the unweighted rms variance and 
then subtract from it the rms variance due to the uncertainties of 
the individual poin ts, which is so-called excess variance (see e.g. 
lEdelson et ai]|2002h ), 

>2 

(5) 
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Note that the variance difference above can be negative if the in- 
trinsic variability is comparable or weaker than the measurement 
uncertainties. If this happens, we set this excess variance to zero. 
We again point out that Sj represent already the fractional rms, i.e., 
it should not be further divided by Gj (which may be zero or nega- 
tive). Then, we calculate the standard deviat ion of the above exces s 
variance, AS], following equation (11) of IVaughan eTail J2003h . 
hereafter V03, 
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We note that the transformation of ASj into ASj is not trivial. V03 
have done it using the standard differential propagation of errors, 
obtaining their equations (B2) and (B3), which, however, we find 
not generally correct. Namely, the assumption behind using deriva- 
tives in propaging errors is that the uncertainty is much lower than 
the estimated quantity. This is often not the case for the excess vari- 
ance, which can be null for either weak intrinsic variability or mea- 
surement errors comparable with that variability, see equation <[5jl, 
whereas its uncertainty is always > 0. Then, the error-propagation 
formula used by V03, ASj = AS] /(dS] /dSj) (using our nota- 
tion), obviously fails, leading to infinite uncertainties. The cause for 
that is the failure of the assumption of ASj <C Sj. To account for 
that, we calculate the uncertainty on the rms without that assump- 
tion, i.e., directly from the defition of the 1-a uncertainty range as 
S]±AS], 



ASj = (S] + AS] 



Si 



(7) 



Here we have chosen the upper error, which is larger than the 
lower one, and which is the only one possible for S] < AS]. 
For AS] < Sj, this becomes the usual ASj ~ AS]/(2Sj) (as 
in V03), which equals ASj ~ aj/K 1 ^ 2 . On the other hand, for 
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Figure 2. Profiles of the orbital modulation of the ASM A (1.5-3 keV, lower crosses) and C (5-12 keV, upper crosses) data at 8 superorbital phase bins. The 
bin number and the phase of the bin center are given on each panel. Note that the minimum and maximum modulation are offset from the and 0.5 phase, and 
appear instead close to 4> = 0.125 and 0.625, respectively. The unit of F is count s . The solid curves give the best-fitting theoretical outflow model (model 
8 in Table[T} described in Section|5] which involves the absorption in the isotropic stellar wind as well as in the bulge situated at the disc edge. The dashed 
curves shows the model component due to the wind only. 



AS? > S?, the result is ASj ~ (ASj ) 1/2 , which should be used 
to correct the upper part of equation (B3) in V03, and which equals 
ASj ~ (2/K) 1 ^ 4 aj (in our notation). Hereafter, we use equation 
10 to estimate the rms uncertainty. 

Note that the above uncertainty estimates are due to the mea- 
surement errors only, and they do not account for the long-term, 
red-noise, variability of the source properties (V03). This is a cor- 
rect procedure for our sample containing most of the currently 
available ASM data, for which we are interested in their actual 
properties, and are not hypothesizing about their behaviour over 
time scales ;§> 10 yr. 



2.3 Hardness ratio 

We also would like to analyse the spectral variability of Cyg X-l 
with orbital and superorbital phase. A useful measure of the spec- 
tral shape is the hardness ratio (HR) of the fluxes in various chan- 
nels, which can be computed in a number of ways. The obvious 
one is to use already available mean fluxes and construct their ra- 
tio. This procedure, however, does not account for short-timescale 
spectral variability. The HR can also be computed for each obser- 
vation (dwell) and then the mean can be obtained. However, we 
have already seen that fluxes follow the log-normal distribution, 
and therefore expect that their ratio could also be distributed in such 
a way. Indeed, Fig.QJb) demonstrates that the logarithm of HR have 
distributions close to normal. Therefore, for the unbiased estima- 
tion of the mean HR and its error, we take the logarithm of HR for 
each observation (using fluxes that are not pre-averaged within the 
P/20 bins) and average them within selected orbital and superor- 
bital phase bins. We also note that the mean HR is computed with- 



out weighting the individual HRs according to their errors, because 
the error is systematically larger for harder spectra (as a result of 
a smaller flux in lower-energy channels), and therefore accounting 
for errors would result in a strongly biased estimate of the mean. 



3 STRENGTH OF THE ORBITAL MODULATION VS. 
THE SUPERORBITAL PHASE 

The folded and averaged profiles of the orbital modulation for the 
ASM A data are shown in Fig. [2] We can see that the orbital mod- 
ulation is variable, e.g., it appears to be the weakest at the super- 
orbital phase $ = 0.625. However, there is also a fair amount of 
statistical noise, and the results of this figure need to be quantized. 
We can see here that the orbital modulation profiles are character- 
ized by rather narrow minima, and thus would not be well fitted 
by a smooth function, e.g., a sinusoid. Thus, we first calculate the 
rms of each dependence to characterize its strength, following the 
method of Section[2] 

Fig-Efa) shows the superorbital phase diagram for the ASM A 
channel. We can see the highly significant flux modulation with the 
superorbital period (cf. L06, Paper I). We also see that the minimum 
of the superorbital cycle is clearly offset from the ephemeris of L06 
by A$ ~ 0.1 (which was based on ~30 yr of data compared to 10 
yr analyzed by us). The crosses in Fig.[3jb) show the corresponding 
rms dependence. We very clearly see a strong dependence of the 
rms on $, with the rms being anticorrelated with the flux. It also 
appears that some phase lag, <0.1, of the maximum of the rms 
with respect to the minimum of the flux is present. 

Fig. 0c) shows the results for all three ASM channels. The or- 
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Figure 3. (a) The superorbital phase diagram for the ASM A channel. The 
unit of F is count s . (b) Comparison of the characterization of the ASM 
A rms dependence using different methods. The crosses are the intrinsic 
rms, S, of the orbital modulation as a function of the superorbital phase. 
The solid histogram gives the amplitude of the orbital variability as fitted 
by sum of three harmonics, see Section [3] The dashed histogram gives the 
corresponding rms for the fitting functions. The solid curves in panels (a) 
and (b) show the dependencies for the theoretical outflow model (model 5 in 
Table[T}. (c) The dependencies of the intrinsic rms of the orbital modulation 
on the superorbital phase for three ASM channels. The crosses with filled 
circles, open triangles and open squares correspond to the channels A, B, 
and C, respectively. 
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Figure 4. (a) The superorbital phase diagram for the Ryle 15 GHz data. The 
unit of F is Jy. (b) The dependence of the intrinsic rms of the 15 GHz orbital 
modulation on the superorbital phase, consistent with being constant. 



larger than the rms simply due to their different definitions. The 
amplitude also shows a strong dependence of the superorbital phase 
similar in shape to that of the rms; however, it appears consistent 
with no phase shift with respect to the flux profile (Fig. [3^). 

We have then searched for a similar effect in the Ryle 15 GHz 
data. We have found, however, that no apparent dependence is seen, 
and the ^-dependent orbital modulation profiles look all similar, 
and consistent with the average orbital modulation (see Fig. 4 in 
L06). Thus, we show here, in Fig. [4] only the results of calculating 
the rms of the orbital modulation as a function of $. In Fig.[4jb), 
we see that the strength of the orbital modulation is consistent with 
being constant, though we cannot rule out some dependence hidden 
in the statistical noise. We have also checked that the 2.25 and 8.30 
GHz data from the Green Bank Interferometer (see L06; Paper I) 
also do not show any statistically significant dependencies. 



bital modulation, due to bound-free absorption, is stron gest in the 
1.5-3 keV range and weakest in the 5-12 keV range i Wen et aTl 
1 19991 ; L06). Consequently, the statistical significance of the depen- 
dence on $ decreases with the energy. 

In order to test the robustness of our finding of the depen- 
dence of the strength of the orbital modulation on $, we have also 
calculated the rms for the ASM A taking into account the weights 
due to uncer tainties of the individ ual points in the orbital phase di- 
agrams (see IZdziarski et aT]|2004l) . This alternative method gives 
only negligible differences with respect to the original one, and 
thus we do not show its results. Then, we have fitted the ASM A 
orbital modulation profiles with a sum of three sinusoidal harmon- 
ics, see equation (2) in Paper I, and calculated both the amplitude, 
(fmax - F,Bin)/ (i^max + F m ia), and the rms for it. In this way, we 
largely avoid contributions to the rms from residual aperiodic vari- 
ability. The results are shown in Fig.[3jb). We see that the values 
of the rms of the fitted functions are very similar to that calculated 
directly from the data in Fig. [2] On the other hand, the amplitude 
(which is sensitive only to the extremes of the fitted function) is 



4 SPECTRAL VARIABILITY 
4.1 Hardness ratio 

The X-ray modulations can also be tracked through the hardness 
ratio. The largest and easily detectable variability is shown by the 
ratio of count rates in ASM channels C and A (C/A). Fig.|5]presents 
the dependence of the mean C/A (computed from the logarithm of 
the ratio, see Section l2~3l l on orbital and superorbital phases. We 
see a strong peak at orbital phase cj> ~ 0, which can be explaine d 
by absorption in the nearly isotropic wind (see I Wen et al .11 19991) . 
The dependence of C/A on the superorbital phase also shows a very 
significant hardening around $ ~ 0. The two-dimensional depen- 
dencies on (j> and $ demonstrate a plateau with C/A~ 1 .4, a signif- 
icant increase in hardness around <j> = 0.0 ± 0.2, and two peaks 
at superorbital phase $ ~ —0.1 and 0.1, which significance is not 
certain. Other hardness ratios C/B and B/A show similar behaviour, 
but of smaller amplitude. 
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Figure 5. The mean hardness ratio in the ASM channels C and A in the hard 
state as a function of the (a) orbital and (b) superorbital phase, (c) Contour 
plot of the smoothed distribution of the hardness ratio C/A over the orbital 
and superorbital phases. 



4.2 X-ray dips 

X-ray dips, which are believed to result from absorption in blobs 
in the stellar win d, are characterised by sig nificant drop in the 
count rate (see e.g . lBalucinska-Church et al. 2000, hereafter BC00; 
iFeng & Cuill2002T) . However, most markedly they manifest them- 
selves by spectral hardening (BC00). It is of interest to study their 
distribution over the orbital and superorbital phase and to compare 
these distribution to the corresponding dependencies of the HR. 

In order to define the dips, we use the ratio of the ASM count 
rates in channels B and A, B/A (HR1 in BC00), and the analogous 
C to B ratio, C/B (HR2 in BC00). We then use the criteria of B/A>2 
or C/B>2.5, which is similar to that of BC00 except that they stated 
that they used both criteria simultaneously. With the present ASM 
calibration, we found only 56 dips satisfying their criterion in the 
ASM data used by us. The cause for almost no dips with both hard- 
nesses large appears to be caused by the dip absorption being par- 




-0.4 -0.2 0.0 0.2 0.4 

Orbital phase, 



Figure 6. The distribution of X-ray dips over (a) orbital and (b) superor- 
bital phase corrected for the coverage. The solid histogram is the hard state 
studied by us, while the dashed histogram is for the entire ASM data set. 
The fraction scale corresponds to the solid histograms, (c) Contour plot of 
the smoothed distribution of all X-ray dips over the orbital and superorbital 
phases. 



tial, i.e., with some small fraction of the flux remaining unabsorbed. 
Then, at a relatively low absorbing column, the flux in the A chan- 
nel is reduced but the B and C channels are only weakly affected, 
so this case yields a B/A ratio increase but not C/B. On the other 
hand, at a column yielding a substantial reduction of the B flux, the 
C/B ratio increases, but the absorbed flux in the A channel is so low 
that it is dominated by the constant unabsorbed component, which 
results in no substantial increase of the B/A hardness. 

In the hard state data, we have found 1151 dips (814 with 
B/A>2 and 387 with C/B>2.5) among 31211 observations, while 
in the whole 10-year data set without any selection we have found 
1336 X-ray dips (995 with B/A>2 and 437 with C/B>2.5) among 
60127 independent observations. Thus most of the dips happen dur- 
ing the hard state. This is expected because the spectral softening 
and increase of the luminosity in the soft state strongly increases 
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the ionization level of the wind, which resul ts in a weaker pho- 
toelectric absorption (BCOO; IWen et al]| 19991) . This also strongly 
confirms the accuracy of our criterion defining the hard state. 

Fig. [6ja) shows the distribution of the dips over the orbital 
phase renormalized to the number of ASM observations in each 
bin. The picture looks relatively similar to fig. 5 in BCOO (based 
on ~2 yr of the data, i.e., five times less than in our data set). The 
peak is at ~ 0, and it is relatively symmetric, especially for the 
hard-state data only. We have also checked that the distributions of 
the dips selected separately in the B/A and C/B look very similar. 
On the other hand, the additional peak at <j> ~ 0.6 claimed by BCOO 
is not found by us, and appears to be due to a statistical fluctuation 
in the previous data set. Indeed, the total number of counts in the 
three bins forming that excess was 34, whereas the continuum level 
(i.e, without the excess) in those three bins corresponds to about 25. 
Thus, the excess corresponds to only ~ 1.5<r in the Poisson statis- 
tics. The existence of the feature at 4> ~ 0.6 is also not supported 
by the dependence of the HR, which shows no signs of spectral 
hardening at this phase (see Fig.[5^,c). 

Then we have studied the distribution of the X-ray dips over 
the superorbital phase. The results are shown in Fig.[6jb). We see a 
maximum around $ ~ 0.05-0.1, which is consistent with the po- 
sition of the flux minimum (see Fig. [3^). The distribution is clearly 
asymmetric relative to the peak, with a slower rise and faster de- 
cline, and it looks like the inver ted flux (i.e., - In F) of Fig.^a). 

Then, the two-dimensional distribution of the dips in <j> and <£> 
is shown in Fig. [6jc). We see that most of the dips that give rise 
to the peak in the orbital phase distribution around <j> ~ 0.0 ± 0.2 
happen around the superorbital phase of $ ~ 0.1±0.2. (The statis- 
tical significance of the presence of two, rather than one, separate 
peaks there is rather low, ~ 2a.) The distributions of the dips re- 
semble strongly that of the HR, which is natural because the dips 
just represent a tail of the HR distribution. 



5 THEORETICAL INTERPRETATION 
5.1 Wind geometry in Cyg X-l 

Let us first summarise our findings. In the radio, we see no mod- 
ulation of the orbital variability with the superorbital phase, while 
in the X-rays such a modulation is visible. In addition to the previ- 
ously known spectral hardening at orbital phase cj> ~ (visible in 
the HR and distribution of the X-ray dips), we find a significant in- 
crease in the HR around superorbital phase $ = 0. This effect can 
be tracked in the dependence of the HR as well as the distribution 
of the X-ray dips. 

Our interpretation of the observed dependencies is as follows. 
The absence of statistically significant superorbital dependence of 
the orbital modulation of the 15-GHz radio emission is consistent 
with the radio being emitted by a jet in the system, in which case 
the orbit al modulation is cause d by wind absorption far away from 
the disc ( Szostek & Zdziarski 2007). For the X-rays, the situation 
is more complicated. 

The X-ray orbital modulation is due to variable absorption by 
the wind of the X-rays emitted close to the disc centre. The absorp- 
tion can be separated into two components. One is independent of 
the superorbital modulation, and is due to absorption in the part of 
the wind steady in the comoving frame, as usually assumed. The 
other component is due to the part of the flow feeding the outer 
edge of the disc, and thus forming a bulge. 

In Cyg X-l system, though the OB star does not fill com- 
pletely the Roche lobe, the wind density is enhanced inside the 



(a) A view along the orbital plane, <t> =0 

observer 




<p=Q 0=0.5 

(b) A view along the orbital plane, * =0.5 

observer 




(c) A view from the top 




Figure 7. A drawing illustrating the effect of a bulge at the outer edge of 
a precessing inclined disc. The material in the bulge absorbs some of the 
X-ray emission originating close to the disc center. The orbital modulation 
due to the bulge is seen to strongly depend on the superorbital phase. In 
addition, there will also be orbital modulation due to the direct wind from 
the supergiant, not shown here for clarity. The elongation of the supergiant, 
almost filling its Roche lobe, is not shown here. A view along the orbital 
plane: (a) the superorbital phase of 0, when the disc is seen closest to edge- 
on and the effect of the bulge is strongest; (b) the opposite case of the su- 
perorbital phase of 0.5. (c) A view from the top, with the arrow showing the 
direction of the observer. The angle (f>^ gives the azimuthal displacement 
of the bulge centre relative to the line connecting the stars, and it is > in 
the case shown here. The maximum of the absorption corresponds then to 
4> = —0b- This view is for any value of <f> except for the shown orientation 
of the elliptical image of the disc, which corresponds to <E> = or 0.5. 



Roche lobe, which is an analog of the Roche lobe overflow 
but by the wind. Suc h a focused wind dFriend & CastoJ 1 19821 : 
iGies & Boltonl [l986b) in some way forms the accretion disc, 
known to exist in the system. The main argument for the existence 
of the disc is an overall similarity of the X-ray spectra and timing 
properties of Cyg X- 1 to those of low-mass X-ray binaries, in which 
case accretion has to form a disc (see, e.g., Zdziarski & Gierlinskil 
l2004h . The disc formation, most likely, leads to a condensation of 
the wind matter near the disc outer edge on the side of the com- 
panion in the form of a bulge, similar to the disc bulge inferred t o 
be present in low-mass X-ray binaries, e. g., White & Holt ( 1982); 
IWhite & Swankl jl999l) ; |Parmar & White! 7i988l) : lHellier& M ason 
ii989|), as illustrated in Fig. [7j On the othe r hand, the bulge can also 
be formed (see, e.g.. lBoroson et al1l200ll) by a shock wave in the 
wind when it encounters the gravity of the companion, the disc, or 
a wind from the disc, which is also likely to be present. In any case, 
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when the fast, > 1000 km s -1 , wind is stopped, the density in- 
creases dramatically. A fraction of the focused wind might pass the 
black hole and be visible as additional absorption at orbital phase 
~ 0.5, however we find no evidence for that in the X-ray data. 



(a) 



observer 



An issue in the above scenario is the position of the bulge rela- 
tive to the line connecting the stars. Consider the accretion process 
in the corotating frame of the binary. In the case of low-mass X-ray 
binaries, the accretion stream leaves the LI point with a small ve- 
locity and, being deflected by the Coriolis force, hits the disc (with 
the outer edge defined by the stream orbital angular momentum) 
at an azimuthal angle 0b ~ 60°, which is measured from the line 
connecting the stars with the origin at the co mpact object (see Fig . 
EJ; and the entry for 4> h - 180° in table 2 in iLubow & Shull 1975b . 
For the mass-rati o in Cyg X-l, q = M BH /M C = 0.36 ± 0.05 
dGies et alj|2003l) , the gas freely falling from LI point would hit 
the disc at 0b ~ 70°. However, these considerations neglect the ra- 
diative acceleration of the stream as well as the diffusive spreading 
of the accretion disc and therefore its potentially much larger size, 
with both effects significantly reducing 0b. 

An additional complexity is brought by a possibility of the 
non-synchroneous rotation of the companion in high-mass sys- 
tems. For example, a slower stellar rotation allows the wind to 
be launched with a non-zero angular momentum in the corotating 
frame and leads to the increase of 0b, while the opposite is true 
for the faster rotation. The ro tation of the compani on in Cyg X-l 
is compatible with corotation dGies & Boltonll986ah . and therefore 
probably does not affect much the gas kinematics. Then, if we mea- 
sure this angle, 0b, in units of the 0-1 orbital phase, absorption of 
the X-ray emission in the bulge will peak at the orbital phase of 
~ 1 — 0b. Indeed, the typical ph ase of major X-ray dips in low- 
mass X-ray binaries is ~ 0.8-0.9 dParmar & Whitel ll988). Some 
other high-mass X-ray binaries show dips at ~ 0.8-0.9, also 
thought to b e caused by the accr etion stream passing through the 
line of sight dBoroson et alj206 H and references therein). 

A crucial further complication in Cyg X-l is that the disc is 
inclined with respect to the binary plane and thus precesses. The 
precession causes changes of the position of the bulge with respect 
to the line of sight. During a single binary revolution, the bulge 
moves up and down, while the inclination of the disc remains ap- 
proximately constant (since P sup ^> P), see Figs. [Tfa, b). At $ 
close to zero, we see the disc at the highest angle, i.e., most edge- 
on. The displacement of the bulge centre 0b relative to the line con- 
necting the stars (see discussion above and Fig.|7J:) will also cause 
a small shift of the superorbital phase at which the bulge absorption 
is maximal. On the other hand, we see the disc close to face-on at 
$ = 0.5, see Fig. |7Jb), and then the bulge is always outside the 
line of sight to the X-ray source. Thus, that additional absorption 
component is absent. 

The above considerations explain the dependence of hardness 
ratio on orbital and superorbital phases as well as the distribution 
of the X-ray dips. Based on the two-dimensional distribution of the 
dips (Fig. |6j;), we have calculated that at least 1/3 of all the X-ray 
dips are caused by the bulge, and the rest are due to the isotropic 
part of the stellar wind. The picture in Fig. [7] can also be used to 
calculate the expected X-ray orbital profiles caused by the wind 
and bulge absorption. We can assume a specific density profile of 
the wind and the bulge, and calculate the optical depth during a 
revolution for a given superorbital phase. 



focused 
wind n • 
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hole 
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observer 
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hole 
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Figure 8. (a) Geometry of the wind, (b) Geometry of the bulge. 



5.2 Model 

Let us consider first the isotropic component of the wind. The wind 
mass density as a function of distance from the center of the star, r, 
can be estimated from the mass conservation law 



Piso(r) 



M 



4:Tvr 2 v(r) ' 



(8) 



where M is the mass loss rate. We assume v(r) oc (1 — i?»/r) ? , 
where i?* is the stellar radius, and consider the attenuation cross- 
section independent of the distance. We thus get the absorption co- 
efficient in the form 



>(r) 



'a 
ao \ - 

> r 



l-R,/c 
l-R*/r 



(9) 



where a is the separation between the black hole and the compan- 
ion, and ao is the absorption coefficient at r — a. We define here 
the characteristic optical depth, Ti SOj o = aao. 

The focused wind can be described by the cone of half- 
opening angle # m ax centred around the line connecting the stars 
(see Fig. [8^ for geometry). The additional opacity can be scaled 
to the opacity of the isotropic component and its angular depen- 
dence can be approxi mated by a parabola dFriend & Castoill 19821 : 
iGies & Bolton|| l986bl) 



a iw (r,8) - a (??fw - 1) 



(10) 



where 8 is the angle measured from the line connecting the stars 
and 7?f w is the ratio of the wind density in the direction of the black 
hole to that of the isotropic component. The total wind absorption 
coefficient is defined by the sum a v (r, 8) = Qi so (r) + ct{ w (r, 8). 
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Let us now compute the optical depth through the wind along 
the line of sight. It depends on the position of the observer. We in- 
troduce the coordinate system centred at the black hole with the z- 
axis along the normal to the orbital plane, and the observer in the x- 
z plane, so that the direction to the observer is n = (sin i, 0, cos i). 
Position of the companion is then a = o(cos <j>, sin <p, 0), where cj> 
is the orbital phase. The angle between the line of sight n and a 
varies with phase: 



cos£ = n • — 



(11) 



The impact parameter is a sin £ and the dist ance of some p oint in 
the wind to the supergiant centre is r = \/s 2 + a 2 sin 2 f , where 
s is its distance to the point of the closest approach (which can be 
negative). The corresponding radius vector is r — n(s + acos£) — 
a (see Fig.[8^). The angle between r and —a is then 



r (-a) 1 , . 2 

cos# = — ■ = - (asm £ 

r a r 



The optical depth through the wind is computed as 

/oo 
a w (r, 9)ds. 
-a cos £ 



(12) 



(13) 



Let us apply this formalism to Cyg X-l. We take the ratio of 
the se paration to the supergiant radius a/R, ~ 2.3 (Ziolkowski 
2005), the inclination i = 40° (see Paper I and references therein), 
the velocity profile exponent ( = 1.05, and the focused wind 
parameters flmax = 20° and r?f w = 3 dFriend & Castor! [T982I ; 
iGies & Bolton] [l986b ) . In this case, the optical depth through the 
isotropic component of the wind from the black hole to the infinity 
in the radial direction away from the companion is about 0.73Ti ao ,o, 
in the perpendicular direction it is 1.26ri ao ,o (i-e. at phases <p — 
0.25, 0.75) and at zero orbital phase along the line of sight r w (0) « 
3ri so ,o. The typical optical depth provided additionally by the fo- 
cused wind across its cone is « ma x(f7fw — l)f Tmo,o » 0.45ri so ,o. 
Thus for Cyg X-l at orbital phases cj> ~ 0, the focused wind adds 
only 15 per cent to the opacity produced by isotropic wind, while 
at <p ~ 0.5, its contribution can reach 60 per cent, but the absorp- 
tion itself at this phase is low. Thus in the first approximation, we 
can use only the isotropic wind model and include the corrections 
introduced by the focused wind later. 

In the case of the bulge, we first need to compute the position 
of the bulge centre, b, relative to the black hole. For the prograde 
precession (see L06) the unit vector along the normal to the pro- 
cessing accretion disc is d = (— sin 8 cos $, — sin 8 sin cos 8), 
where 8 is the precession angle. Assume now that the bulge centre 
lies at the disc plane and the projection of b on the orbital plane 
x-y makes an angle (p b with the line connecting the black hole to 
the companion (i.e. the azimuth of b is <j> + (p b , see Fig. [7};). We 
then get the unit vector of the bulge centre 

[cos(<ji + <j} h ),sin(<j) + 4> b ),ta,n5cos((j> + <j) h - $)] 

O = . (14) 

y/l + tan 2 <5cos 2 (^ + <p b - W) 

The angle it makes to the line of sight is given by (see Fig.[8j)) 

cos = b ■ n (15) 
sini cos(4> + <p b ) + cos i tan 8 cos(0 + <p b — <t) 
s/1 + tan 2 <5 cos 2 (0 + <j) b - $) 

Let us assume an exponential dependence of the absorption 
coefficient on the distance p from the bulge centre, 



with rb being the bulge scale-height. This gives the optical depth 
from the bulge centre to infinity of Tb,o = rbib,o- On the other 
hand, the optical depth from the black hole through the bulge along 
the line of sight, 



a b (p)ds, 



(17) 



iicos/3 



depends on the orbital as well as superorbital phase. Here R is the 
distance to the bulge centre from the black hole (i.e., approximately 
the disc size) and p = \J s 2 + R? sin 2 j3. 

For simplicity we assume that the wind and the bulge are in- 
dependent and therefore the orbital modulation profile is given by 



F(ct>,$) = F (cosip) exp[-T„((j))} exp[-r b (^, $)], 



(18) 



where Fo is the intrinsic flux (which depends on <J>) without ab- 
sorption in the direction of the observer and ip is the angle between 
the disc normal and the line of sight: 



cos ip = n ■ d — cos i cos 8 — sin i sin 8 cos $. 



(19) 



a b (p) = a b ,oexp(-p/r b ), 



(16) 



The retrograde precession can be modelled by substituting $ 
— $ in the above formulae. 



5.3 Modelling the data 

In order to describe the profiles presented in Fig. [2] with the model 
of Section |5T2l we need to specify the angular distribution of the 
intrinsic flux, Fo (cos ip) . In Paper I we have considered four simple 
analytical models: 

(a) the black body, with the flux proportional to the projected 
area, Fo (cos tjj) = A cos ip. 

(b) an anisotropic model of Fo (cos ip) = Acostp(l + 
rjcostp) with parameter r\ giving the degree of deviation from 
the black body. Such anisotropy ca n be produced for example 
by thermal Comptonization (Paper I; ISunvaev & Titarchuk|[l985l 
IViironen & PoutanerJ 1 20041) . which the dominant radiative pro- 
cess giving rise to X-rays in the hard state of Cyg X-l (e .g., 
iGierlinski et al.lll997l : |Poutanen & Coppilll998l : Foutanenlll998l) . 

(c) the steady jet model, Fo(cosi/>) = -A[7j(l — 
Pj cos ip)] - ( 1+r \ where /3j = v/c is the jet velocity, 7j = 

1 / yl — P 2 is the jet Lorentz factor, and Y is the photon in- 
dex of the X-ray radiation. By the 'jet', we mean here either the 
base of the jet in the dire ct vicinity of the b l ack hole, or an out- 
flowing corona (se e e.g. iBeloborodovl 1 19991 : iMalzac et al 1 l200ll : 
iMarkoff et a"i]|2005h . 

(d) the slab absorption model, Fo(cosip) = 
A cxp(— t/ cos ip), which can be associated, for example, 
with some kind of a disc outflow. 

All the models provide a good fit to the superorbital variability 
of Cyg X-l (Paper I). Models (b) and (c) can be considered as more 
physically motivated, but we consider here all of them. In order to 
keep the number of parameters to minimum we fix the inclination 
of the system i = 40°. The precession angle 8 is not well deter- 
mined in models (b)-(d) as it is anticorrelated with other parame- 
ters {q, /3j , r, see Paper I). Thus we fix it at three values between 5° 
and 10°. 

The parameters describing the absorption of radiation are the 
characteristic optical depths Ti so ,o and r b ,o for the wind and bulge, 
respectively. Additional parameters are the bulge density scale 
measured in units of the disc size, r b /R, and the phase shift, <p b , 
of the position of the bulge centre. An arbitrary shift in the su- 
perorbital phase, A$ (due to the uncertainty of the superorbital 
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Table 1: Best-fitting model parameters. 
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ephemeris), is also introduced (i.e. we replace $ by $ — in 
all formulae of Section \52\ . The parameters describing the radia- 
tion pattern are the normalization, A (for ASM A channel), and the 
anisotropy parameter, r), in model (b), /3j in model (c) (where we 
fix T at a typical hard-state value of 1.7), and the slab optical depth, 
r, in model (d). 

We consider the prograde precession of the disk (L06). In or- 
der to understand the influence of the model complexity on the re- 
sults, we consider first only the isotropic component of the wind 
(models W in TableQ} and fit the ASM A profiles only, which show 
strongest variability. We find that parameters r^/R and Tb,o are 
anticorrelated, and cannot be determined separately. This happens 
because various combinations of the two parameters can give the 
same optical depth through the bulge at a given impact parameter. 
Therefore, we fix rb/R = 0.2. The best-fitting model parameters 
are presented in TableQ] For model (a), the precession angle agrees 
within the errors with the results of Paper I. The jet model, (c), pro- 
vides a slightly better fit for smaller precession angles. The mod- 
els (b) and (d) also give statistically similar fits. The phase shifts 

and 0b are well constrained by all the models. The fits require 
the shift of the bulge centre from the line connecting the stars by 
4>b ~ 0.07 (i.e., 25°, see Fig. [7};). All the models give similar opti- 
cal depths through the wind and the bulge. The wind optical depth 
r w varies between 0.28 (i.e. ~ 3ri SOj o) and 0.08 (i.e. ~ Ti SO; o) for 
tj> varying between and 0.5. For the bulge, Tb varies between 0. 15 
and 0.007 at $ = and between 0.05 and 0.008 at $ = 0.5. 

We now add an additional focused wind component with the 
parameters specified in Section |5T2l and fit the data using jet model 
(c). The resulting best-fitting parameters are not very much differ- 
ent from those obtained with the isotropic wind model (compare 
entries 5 and 7 in Table 1). This is expected, because the focused 
wind affects the total opacity on average at about a 30 per cent level. 

Finally, we fit the light curves in channels A and C simultane- 
ously. Two additional parameters have to be introduced: the ratio of 
the absorption coefficients (and optical depths) in channels C and 
A, tc/ta, and the ratio of the normalizations (intrinsic hardness 
ratio), C/A. The best-fitting results for the main model parameters 
change only slightly (compare entries 5 and 8 in Table 1). Because 
the mean absorption coefficients in channels A and C differ only by 
a factor of 3, the absorbing gas has to be rather strongly ionized. 

For the retrograde precession, all these models give much 
worse fits to the data. 



6 DISCUSSION 

6.1 The origin of beat frequencies 

A collateral effect of the coupling between the orbital and super- 
orbital modulations may be appearance of additional frequencies 
in the power spectrum. If the two modulations were independent, 
there would be simply two peaks in the power spectrum at the cor- 
responding frequencies. On the other hand, if one modulation de- 
pends on the other, beat frequencies, at v — 1/P ± 1/P sup , may 
appear. Indeed, L06 reported finding the lower of the beat frequen- 
cies (albeit at a relatively limited statistical significance), and also 
found that its origin from X-ray reflection from the surface of the 
companion is unlikely. 

Here, we have tested whether the discovered dependence of 
the orbital modulation on the superorbital phase may indeed cause 
beat frequencies to appear. Using our model (given by equation 
dl8t and other formulae of Sections 15.21 15.31 with parameters of 
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Figure 9. Power density spectra predicted by our models (arbitrary nor- 
malization). The solid curves show the power spectrum of the flux for the 
outflow model, 5 in Table [T] The inset zooms on the frequency range near 
l/P. The dashed curves show the model with the absorption only in the 
wind, i.e., neglecting the presence of the bulge. 

model 5 in Table [TJ, we have generated a light curve and com- 
puted Fourier power-density spectrum (PDS). We have found that 
our model gives rise to strong peaks at frequencies 1 / P Bup and 1 / P 
with harmonics as well as to two peaks in the power spectrum at the 
beat frequencies. Interestingly, the lower beat-frequency peak is 3.7 
to 5.4 times stronger than the higher one (depending on whether we 
compute Fourier transforms from the flux or from the logarithm of 
the flux). Fig. [9] shows the flux Fourier transform for this case for 
the outflow model. We then compare these predictions with a sim- 
pler model where absorption in the bulge is neglected. In this case, 
there are two beat-frequency peaks of equal strength in the PDS of 
the flux (see the dashed curves in Fig.[9j>, while they are missing in 
the PDS computed from the logarithm of the flux, because the cou- 
pling disappears. If on the other hand, only bulge produces absorp- 
tion, PDS shows both beat-frequency peaks with the strength ratio 
of 10 and 5 for the flux and its logarithm, respectively. Finally, we 
experimented with the model where intrinsic flux, Fo, as a function 
of superorbital phase was assumed to be constant and both bulge 
and wind are responsible for absorption. Now, the strength of the 
peak at 1/ P sup has diminished by three orders of magnitude, while 
the behaviour of PDS at l/P and the beat frequencies was almost 
identical to the full model with variations of Fo (the ratio of peak 
strengths is 6.5 and 5.5 for the flux and its logarithm, respectively). 

We see that the lower beat-frequency peak is always stronger 
than the higher one when absorption is modulated by the bulge 
(for prograde precession). The coupling discovered in this work 
thus predicts a presence of beat frequencies with a stronger low- 
frequency peak, which is consistent with the discovery by L06 of 
only the low-frequency peak. 

6.2 Superorbital variability and outbursts of Cyg X-l 

It is of interest to consider whether the superorbital variability of 
Cyg X-l is related to other aspects of the sourc e activity. Recently, 
the MAGIC collaboration d Albert et at] 1 20071) reported detecting 
TeV emission from Cyg X-l. That detection, on MJD 54002, 
took place in the m iddle of a strong X-ray outburst of Cyg X-l 
( lTurleretalJl2006h . We have checked that that time corresponds 



to the peak of the superorbital cycle, $ ~ 0.5, when the disc 
and jet of Cyg X-l are most face-on. On the other hand, L06 
found that the superorbital cycle was uncorrelated with the ap- 
pearance of other strong X-ray outbursts of the s ource of duration 
of da ys dStern. B eloborodov & Pouta ner] 1200 ll ; iGolenetskii et al.l 
|2003j). Thus, the significance of the coincidence of the TeV burst 
with the peak of the superorbital cycle in the present case remains 
unknown. 

Inter estingly, the orbit al phase of the TeV outburst was at 
4> ~ 0.9 ^Albert et al.|[2007l) , at which absorption of TeV photons 
by pair production on the stellar photons is very strong. A possible 
way to obtain detectable TeV emission is then via pair cascades. 
We note that the statistical significance of the detection was rela- 
tively limited, 4.1<r, and thus an independent confirmation of this 
deteciton is desirable. 



7 CONCLUSIONS 

We have discovered the dependence of the orbital modulation 
strength and the hardness ratio on the superorbital phase of Cyg 
X- 1 . The observed effects can be explained by the presence of the 
absorbing material more or less fixed in the corotating frame of 
the stars. We associate this material with the bulge formed by the 
accreting stream impacting the accretion disc. Because of the disc 
precession (causing superorbital variability), the bulge moves up 
and down and its influence on absorption varies. At the superor- 
bital phase 0.5, the line of sight does not pass through the bulge, 
while at $ ~ the absorption in the bulge is maximal. We es- 
timate the maximal optical depth at 1.5-3 keV through the bulge 
(for our line of sight) of about 0. 15, while the stellar wind produces 
twice as much of the absorption. 

Using a simple model of the bulge and the stellar wind in- 
corporating the angular dependence of the intrinsic X-ray radiation 
from the black hole vicinity, we were able to reproduce the detailed 
shape of superorbital variability as well as of the orbital modulation 
at various superorbital phases. We find the bulge centre is displaced 
from the line connecting the stars by about 25° . We also study the 
distribution of the X-ray dips over superobrital phase we find their 
concentration towards the superorbital phase 0.1, which coincides 
with the position of the flux minimum. We thus are in position to 
claim that the X-ray dips observed in Cyg X-l at around zero or- 
bital phase have direct relation to the bulge which, in turn, causes 
variation of the orbital modulation with the superorbital phase. We 
Fourier analyse our model, and find it explains the finding of only 
the lower beat frequency between the orbital and superorbital fre- 
quencies in the observed power spectrum (L06), provided the disc 
precession is prograde. 

We also find that both the X-ray and radio fluxes of Cyg X-l in 
the hard state on time scales >10 4 -s have lognormal distributions, 
which complements the finding of a lognormal flux dis tribution in 
the hard state on ~l-s time scales dUttlev et al.| [2005). We stress 
out that the lognormal character of the flux distribution requires 
that flux logarithms rather than fluxes themselves should be used 
for averaging and error analysis. We also correct a mistake in the 
treatment of V03 of the uncertainty of intrinsic rms variability of 
light curves in the case when the uncertainty is higher than the in- 
trinsic rms (which is often close to null). The mistake stems from 
the failure of the assumption of the uncertainty to be much less than 
the estimated quantity, used in the standard propagation of errors. 
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